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1 Introduction 



The appearance of a light Higgs boson, which might be similar to that of the Standard 
Model (SM), is a specific feature of the minimal supersymmetric Standard Model (MSSM), 
and distinguishing the various scenarios is an important task for coming experiments at 
hadron colliders. Higgs discovery through the main decay channel into a pair of bb quarks 
will be difficult [1], in particular in the mass range below 140 GeV which is common to 
both the SM and the MSSM, and experimental investigations have to deal with signatures 
resulting from the more rare decay modes. 

Besides inclusive single Higgs production, Higgs boson production in association with a 
high-pT hadronic jet provides a useful channel for Higgs searches at hadron colliders, which 
allows for refined cuts increasing the signal-to-background ratio. Specifically, for the SM 
Higgs boson [2], simulations considering the decay channels ^ 77 [3, 4] and H — > t~^t~ 
[5] have shown promising signal and background results for the ATLAS detector. 

In the meantime, a lot of progress has been made towards improving the SM pre- 
dictions. The fully differential distribution for Higgs production at next-to-next-to-leading 
order QCD accuracy has become available [6], improved by resummation of logarithmically 
enhanced terms for low px [7] . Higher-order corrections to differential cross sections for a 
Higgs boson associated with a high-p^ jet have been obtained explicitly: the next-to-leading 
order QCD corrections in the large top-mass limit [8] and, recently, the corresponding re- 
summation of soft-gluon emission effects [9]. 

For the analogous MSSM processes, pT-distributions have been studied in the limit of 
vanishing superpartner contributions at leading order [10] and were improved recently by 
soft-gluon resummation effects [11]. This limit is usually a good approximation when the 
superpartners are heavy, at a mass scale around 1 TeV. The general case was treated in a 
previous paper, where we studied the MSSM prediction for the total hadronic cross section 
for the production of the lightest neutral Higgs boson /i° in association with a high-pr 
jet including all superpartner loop contributions at leading order [12], thus extending the 
MSSM predictions also to the case of light supersymmetric partners of quarks and gluons. 

In the present paper, we complete our study by presenting results for the distributions 
of pseudorapidity and transverse momentum, 773 and pr- Section 2 briefiy reviews the 
contributing parton processes, and Section 3 introduces the expressions for the differential 
hadronic cross sections da/drj^, da/dpT and d^a /drj^dpT in the presence of cuts. In Section 4 
we present numerical results for the distributions at the LHC and the Tevatron, based on 
MSSM parameters around specific benchmark scenarios, in comparison with the SM results. 

2 Partonic processes 

At the partonic level, production of a Higgs-boson hP together with a jet is described 
generically by 

Pi{k^)+P2{k2)^PM+h%k^), 
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involving partons Pi with momenta ki — {Ei,ki). There are three classes of parton pro- 
cesses: gluon fusion g + g g + (Figure 1), quark-gluon scattering q + g ^ q + hP 
(Figure 2), and quark-ant i- quark annihilation q + q ^ g + (Figure 3). While gluon 
fusion is an entirely loop-induced process, the other two classes also get contributions from 
tree-level 6-quark initiated processes (see Figure 4). Those Born- type processes are in gen- 
eral dominant for 120 GeV, while for large values of uia the loop-induced processes 
dominate. This behaviour is essentially a consequence of the Yukawa-coupling of the light- 
est MSSM Higgs boson to 6-quarks, which can be enhanced for low values of the A-boson 
mass TTiA- Formulae for the relevant couplings occurring in this process in our notation can 
be found in [12, 13]. 

The virtual presence of the superpartners in the loop contributions modify in addition 
the overall production rates for supersymmetric Higgs bosons. Moreover, they affect the 
angular distributions and thus, at the level of hadronic processes, change rapidity and 
transverse-momentum distributions of the Higgs bosons or the jets, respectively. 

In terms of the invariant kinematic variables 

■s=ik^ + k,f, i=ik,-ks)\ (1) 

the partonic differential cross section reads 

da^^ 1 



with the spin- and colour-summed/averaged ^-matrix element squared. 

Instead of t, the cross section can be expressed in terms of other variables, such as the 
transverse momentum pr — k^^T, or the scattering angle 9 or rapidity ys in the parton 
center-of-mass system (cms), yielding the representations 

da ,^ , da^ ,^ , da^,^ , 

-i—\^^Pt) = -, — {s,pt) + - — (s,pr), (3) 
apT dpT dpT 



da^^^\^ , 2\ki\pT da,^ 

— {s,pt) = . :^ 



t=tF[S](s,PT) 



do' ^1^11^1. da , ^ 

— (s,^) = 2|A;i||A;3|sm^ ^{s,t) 

d9 dt 

— s,7/3 =2 A;i E3 l-tanh^T/a s,^! 
dy-i dt 



t=t{s,e) 



with 



(5) 
(6) 



tF[B]{s,PT) = -2E1E3+ H 2\k,\^ \k3\^ - pS , (7) 

i{s,m) = -2^1 ^3 + 2|fci|^3tanh^3, (8) 

i{s,9) = -2£;i£;3 + 2|A^i||4|cos^, (9) 

E, = \k,\ = 4- £3 = |£,i = if#. (10) 
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3 Hadronic observables 



Hadron colliders like the LHC (Tevatron) collide protons with protons (anti-protons) in 
the center-of-mass system (the laboratory frame) with a total energy and individual 

— * — * 

momenta P and —P, respectively. Hadronic cross section are obtained via convolution 
of the parton-level cross sections with the parton distributions and summation over the 
various contributing partons. Experimental restrictions to the detectability of the produced 
particles arc conventionally realized by imposing specific cuts to the kincmatically allowed 
phase space. Typically, cuts are imposed on the final-state transverse momentum and/or 
the pseudorapidity in order to have high-p^ jets not too close to the beam axis. In our 
case, we choose the following selection criteria, 

Pt = \h\ sm9 > p^'" , |?73| < ?7max , (11) 

where pt and 773 denote transverse momentum and pseudorapidity of the final state parton. 
The first condition leads to an energy cut on the invariant mass of the parton system, 



VI > pp'' + J ml + {pf-y = ^/I'o , (12) 



and to a cut on the scattering angle 9 in the parton cms. 



9 A „min 

sin^>f^. (13) 



ml 



The second condition in (11) avoids partons with angles 6\sh unobservably close to the beam 
axis in the laboratory frame. For an outgoing massless gluon or quark with momentum k^, 
773 coincides with the parton rapidity ^3^'' in the laboratory frame, which is related to the 
rapidity in the parton cms via — y^^ + iplr, x), where ip is determined by the relative 
velocity 

/3= (x^ -r)/(x2 + T) = -tanhV'(T,x) (14) 

of the parton cms with respect to the laboratory frame. Thereby, it is assumed that the 
initial state partons carry momentum fractions x and t/x oi the hadrons they originate 
from. The condition I773I < 77max can be written as follows, 

arccos(tanhr7max) < ^lab < tt — arccos(tanhr7max) ■ (15) 
This translates into the following additional cut on the scattering angle in the parton cms, 

^min < 6* < ^max (16) 

with 

^^^^ _ (r - x"^) + (r + x"^) tanhr^max 

(r + a;2) + (r - a;2) tanhr^max ' 
cos^ ^ ~ ~ + tanh77in-^ 
™^ (r + x2) - (r - a;2) tanh T^^ax ' 
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Thus, in the additional presence of the pseudorapidity cut the partonic cross section 
for initial state partons n and m, with momenta xP and —{t/x)P, depends on rjjnax 
and both variables r and x, 



(19) 

s=tS 



3.1 Integrated cross section 

The hadronic cross section for hP + jet production from two coUiding hadrons, A and B, 
contains the parton densities and the partonic cross sections anm from (19) as follows [14], 



(7AB('S',77max,PT"') = / / dx ] — X 



) X 

X 

+ fm/Aix, IJ'F)fn/Bi-, A«f) x,Vm^,pD}, (20) 

X ) 

where fn/A{x, fJ'F) denotes the density of partons of type n in the hadron A carrying a 
fraction x of the hadron momentum at the factorization scale fip [{A,B) = {p,p) for the 
LHC and {p,p) for the Tevatron]. The sum over unordered pairs of incoming partons 
runs over {n, m} = gg, qg, qq with the outgoing parton g, q, g respectively, for the various 
channels. The lower bound of the r-integration (tq) is determined by the minimal invariant 
mass of the parton system, sq = ^oS, according to (12). 

Care has to be taken for the fact that a forward-scattered parton n out of hadron A 
combines with a backward-scattered parton n out of hadron B or vice- versa to a given value 
of ^lab or 773. Hence, a distinction has to be made in the notation for the integrated partonic 
cross section: anm has the parton n moving in the direction of the incoming hadron A in 
the laboratory frame and amn has it in the opposite direction. 

3.2 Transverse- momentum distribution 

The hadronic cross section differential in pr is given by the convolution integral of the 
corresponding partonic differential cross section^. 



dpT f Jro Jr X [ l + Onm 

{n,m} " 

■{tS.pt) e(r?max - Kf(t,2;,Pt)|) 



X 



{n,m} 

d&l 
dpT 



B 



+ \ m^n\ , (21) 



^In the following we drop the factorization scale iip the notation. 
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with 

(r, ^, Pt) = -Hr, x) + [-]arcosh ^^^^ , %7 = ,73"- , 7^3-^ = %7 , 
and -0 from (14). 

3.3 Pseudorapidity distribution 

Taking into account the two directions of motion of the initial state partons in the labora- 
tory system it is convenient to define 

vTim, T, x) = + V'(t, x) , yr(V3, r, x) = -yr(V3, r, x) . (22) 

With this convention the formula for the hadronic cross section differential in 773 reads as 
follows, 



{n,m} "^0 T 
fn/A{x)fm/B{l) da 

1 



1 + 5nm dys 

where 



i^S, j/3) 



e(pT(T, X, 7^3) - pT"") + (n^m)} (23) 



Pt{t, X, 7]3) = — ^ . (24) 

2\/tS cosh(773 + -0(2;, r)) 
3.4 Two-fold differential cross section 

The hadronic cross section differential in pt and 7/3 can be written as an integral over 
a single parameter, which can be chosen to be the parton momentum fraction x. The 
expression follows from the representation 

^{S^m.PT)- (25) 

dpTdr]-i 

^ dr / — Q{x -T)Gnm{x,T,r]3) 5{pT{T,x,r)3) -Pt)) , 

with 



^ / V fn/A{x)fm/B{-) da nm , r< - ^ 

Gnm{x,r,ri3) = — --p—{TS,ys) 

1 + dnm dys 



+ {n^m), (26) 

^3=^3™-('?3,T,a;) 
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by performing the r-integration with the help of the (5-function, yielding 

(faAB 



E 

{n,m} 



(27) 



^dx 

X 



e{x- T)e(T - ro)^=- y= —Gnm{X,T, 7]^) 



T=Ti(x,ri3,pT) 



where 



Ti{x,'r]s,pT) = X 



xpre + ^ 
x\fS — pt e^'-> 



(28) 



fulfils the relation 



PT{x,Ti,r]3) -pT^O 



(29) 



with Pt{x, t, rjs) from (24). 



4 Numerical results 

In the following discussion we want to illustrate the MSSM predictions for the pseudorapi- 
dity and transverse momentum distributions of the hadronic processes pp — > + jet + X 
and pp ^ + jet + X and outline differences between MSSM and SM predictions. For 
comparison of a given MSSM scenario with the SM we choose the SM Higgs mass to have 
the same value as the /i" boson in that MSSM scenario. For the numerical evaluation we 
use the cuts (11) with = 30GeV and r^max = 4.5 as standard cuts, which have been 
used in previous Standard Model studies for the LHC [3, 4]. 

The evaluation has been carried out with the MRST parton distribution functions [15], 
with the renormalisation scale fiji and factorisation scale for the gluon and the light quarks 
IJ,p^ chosen both equal to mh. For the bottom-quark factorisation scale /^^^ we choose 
nih/A, inspired by the NNLO prediction for the process bb — > hP [16] where it proved to 
be the proper scale choice, as anticipated by several authors [17]. For the strong coupling 
constant 0:5 (///?), we use the formula including the two- loop QCD corrections (see e.g. 
Ref. [18]) for n/ = 5 with A^^^ = 204.8 MeV. 



4.1 Parameters 

We adopt for our discussion the MSSM benchmark scenarios for the Higgs search at LEP 
suggested in [19] except for the large-// scenario which has been ruled out by LEP data [20]. 
We are interested in effects from the virtual superpartners. Therefore, the no-mixing and 
j^max scenario are generalised to have a lower common sfermion mass scale MgusY than in 
the original proposal. The two remaining scenarios in [19], the gluophobic and small-aeff 
scenario, have already a rather low sfermion mass scale. 
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Although our aim here is mainly exemplary we try to take into account relevant pa- 
rameter constraints from previous experiments. 

Firstly, we calculate for each parameter point the MSSM predictions for m^o and 
a{e'^e~ — >■ h^Z) x BR(/i° — > bb) and exclude it if the m/jo- dependent LEP-bound on a x BR 
is violated (according to Table 14(b) of [20]). We use FeynHiggs 2.5.1 [21] for calculating 
the m/jo prediction and allow for a theoretical uncertainty of 3 GeV. 

Secondly, wc calculate the leading order MSSM prediction for the branching ratio 
BK{B — >■ Xg^f) [22] and exclude parameter points if the prediction falls outside of the 
range (3.55 ± 1.71) ■ 10~^. This range is determined by using the experimental central 
value [23] and adding up the experimental 3a interval {p^ 10~^) and an estimate of the 
independent theoretical uncertainty (0.71 • 10~^). The latter estimate (20%) is guided by 
the detailed discussion of theoretical uncertainties for the SM prediction [24]. 

Furthermore, we took care that all mass exclusion limits from direct search results 
for superpartner particles [18] are met and that the dominant squark-contribution to the 
electroweak p-parameter[25] stay within ±0.0025. 

Interestingly, applying the rather conservative bound derived from the BR{B — > Xg^) 
prediction, it turns out that the gluophobic scenario is ruled out in the range we study, 
{ruA e [50 GeV, 1000 GcV], tan/? G [1, 50]}. For this scenario BR{B X,-f) is notoriously 
too large with typical values of the order of 10~^ and the lowest value being 7 • 10"''. We 
do not include results for this scenario here. 

A brief specification of the scenarios investigated is following. 

no-mixing{700) scenario : The soft-breaking sfermion mass parameter is set to Msusy = 
700 GeV. For Msusy significantly below 700 GeV the whole range of the m^-tan /3- 
plane which we study is ruled out. The off-diagonal term Xf {— At — //cot/3) in 
the top-squark mass matrix is zero, corresponding to a local minimum of rrih as a 
function of Xf. The supersymmetric Higgsino mass parameter fi is set to —200 GcV, 
the gaugino mass parameters to Mi = M2 = 200 GeV, and the gluino mass to 
Mg = 800 GeV. When tan P is changed. At is changed accordingly to insure Xt — 0. 
The settings of the other soft-breaking scalar-quark Higgs couphngs are A^ — At and 
Aq — {q — u, d, c, s). 

The input for the Higgs sector is specified by uia — 500 GeV and tan /3 = 35. 

^max^^QQ^ gQQjiQj-lQ ■ Xt is set to 2MsusY which yields the maximal value of mh with 
respect to stop mixing effects. We set Msusy = 400 GeV and the other parameters 
are chosen as in the previous scenario. 

In the m5^^(400) scenario small values of niA are still allowed. Hence we examine two 
Higgs sector scenarios: rriA = 110 GeV, tan/? = 30, and rriA = 400 GeV, tan/? = 30. 
The former leads to the dominance of 6-quark initiated processes, while the latter is 
dominated by the loop-induced processes [12]. 

small-aejf scenario : This scenario gives rise to suppressed branching ratios for the decays 
— > bb and t+t~, especially for large tan /? and moderate values of m^. The settings 
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are: MgusY = 800 GeV, Xt = -1100 GeV, Mi = M2 = 500 GeV, // = 2000 GeV and 
M~g = 500 GeV. 

We choose = 400 GeV and tan = 30 in the Higgs sector. Values niA ^ 300 are 
ahcady ruled out, and therefore all scenarios with dominance of the 6-quark initiated 
processes as well. 

4.2 Differential cross sections at the LHC 

The crucial parameter determining the properties of /i°+jet production in the MSSM is 
tua [12]. For niA ^ 120 GeV and tan/5 not too small (> 5) the 6-quark initiated processes 
(see Figure 4) dominate the production rate by far, while for larger ttia this role is taken 
over by the loop- induced processes (see Figs. 1 to 3). Accordingly, we split our discussion 
into the high-m^ and low-m^ cases. 

4.2.1 High rriA 

The influence of rather light, yet not excluded, supcrpartners on the total hadronic cross 
section has been demonstrated to be strong [12]. In particular for the m™'^^(400) scenario 
with MsusY = 400 GeV the MSSM cross section for tua > 200 GeV and any tan /3 G [1, 50] 
is reduced by about 20 — 40% compared to the SM. Here, we investigate the impact on the 
shape of the differential distributions with respect to the SM. 

No-mixing(700) scencirio: 

Figure 5 displays the results for the no-mixing (700) scenario. The lower left panel shows 
da/dris for the MSSM and SM process and also for the three types of subprocesses con- 
tributing to h^+iet production individually. The typical two-peak shape of the quark-gluon 
scattering contribution is caused by the harder momentum distribution of up- and down- 
type quarks with respect to their anti-particles, which lead to a net boost in the direction 
of motion of the proton providing the quark, still visible in the sum over all partons. The 
gluon- fusion and the small qq contribution are peaked around rj^ = 0. 

The upper left panel of Figure 5 shows the relative difference, 5, between the MSSM 
and SM prediction for da/drj^. While the total hadronic cross section in the MSSM is 
enhanced by about 6% compared to the SM, the enhancement of da/drjs varies only by 
about 1.3% (between 5.2% and 6.5%) in the range |?73| < 4.5. 

The right panels in Figure 5 show da/dpT for the MSSM and SM processes and the 
corresponding relative difference, 6. The thickness of the lines in the lower right panel 
hides the few-percent deviations between MSSM and SM for the important contributions 
from gluon-fusion and quark-gluon scattering. Interestingly, the deviation between MSSM 
and SM is largest for qq scattering, e.g. 14% for — 50 GeV. The relative difference in 
da/dpT between the MSSM and the SM (Figure 5, upper right panel) varies between 3% 
for pT = 50 GeV and 8% (17%) for pr = 500 GeV (1 TeV) 

The relative difference between the MSSM and SM prediction for the two-fold dif- 
ferential cross section d'^a/dpT/drf^, indicated by the contours in Figure 10(a), shows a 



8 



quantity 


no- mixing (700) 
SM MSSM 


SM MSSM 


small aeflt 
SM MSSM 




1.623 pb 1.762 pb 
1.682 pb 1.749 pb 


1.448 pb 1.096 pb 
1.419 pb 1.031 pb 


1.490 pb 1.356 pb 
1.480 pb 1.299 pb 


R = ac/af 


0.965 1.008 


1.020 1.063 


1.007 1.044 


A 


4.5% 


4.2% 


3.7% 



Table 1: Cross section predictions for Higgs + jet production with jets radiated into the 
central ((Jc) and forward part of the detector (cx/), together with their ratio R and the 
relative difference between the MSSM and SM value for R, A. 

non-trivial behaviour with an overall variation of more than 6% in the depicted range, 
I773I < 4.5 and 30GeV < pt < 500 GeV. The differently shaped dots in Figure 10(a) show 
the absolute difference between MSSM and SM, which gives an idea of the kinematical 
region where the LHC experiments may become sensitive to this difference. 

Modifying the cuts may increase the sensitivity to deviations from the SM. Guided by 
Figure 10, we study the cross section cr/ with rather soft forwardish jets and cTc with harder 
more central jets: 

(7c = cr {pp h^+jet + X) ||r,|<1.5,pT>70GeV , (30) 
CT/ = cr (pp /l° jet + X) |i.5<|r,|<4.5, 30 GeV<pT<50 GeV ■ (31) 

The results are put together in Table 1, where also the ratio 

R^^ 

and the relative difference between MSSM and SM 

^ _ -Rmssm — RsM 
RsM 

are listed. While each individual cross section in the MSSM and the SM is still of the order 
of 1 pb, which translates into 10^ raw events for an integrated luminosity of 100 fb~^, the 
MSSM ratio -Rmssm differs by 4.5% compared to Rsu- 
m™'*''(400) scenario: 

Figure 6 contains da/drjs and da/dpr for the m™^(400) scenario and the corresponding 
relative differences to the SM prediction. The total hadronic cross section in this MSSM 
scenario is about 25% smaller than in the SM. Yet, as far as the 773 and pt dependent 
differences between MSSM and SM are concerned, the same qualitative picture appears. 
The variation of the relative difference 6 with rj^ in the range |?73| < 4.5 is about 2% and 
with Pt in the range pr e [30 GeV, 1000 GeV] is about 7%. 

Similar to the no-mixing(700) scenario, the difference in the doubly differential cross 
section in Figure 10(b) has a non-trivial rjs and pr dependence. This suggests a similar 
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(32) 



(33) 



refinement as in the no-mixing(700) scenario. For simpficity, we calculate the same ratio 

of cross sections according to the same cuts as before (see Eqs. (30) to (32)) for the SM 
and MSSM. Table 1 shows A — 4.2% as the relative difference of the ratios. 

Small-CKeff scenario: 

In Figure 7, da/drj^ and da/dpT for the small-aeff scenario and corresponding relative 
differences to the SM are shown. While the total hadronic cross section in this scenario is 
about 11% below the SM, the variation in the 773 distribution is about 2%, as in the other 
two scenarios. Opposite to the other two scenarios the pr spectrum is slightly softer than 
in the SM. The range of variation of da/dpT with p^ is about 15%. 

The doubly differential cross section in Figure 10(c) shows a behaviour similar to the 
other two scenarios. Calculating the cross section ratio R according to Eqs. (30) to (32) 
we get A = 3.7% (see Table 1). 

4.2.2 Low tua 

As an example for the low-m^ case at the LHC we show results for the m™^(400) scenario 
in Figure 8. The change with respect to the SM is dramatic. Due to the enhanced cross 
sections of the 6-quark processes, the quark-gluon scattering contribution dominates the 
cross section and even the contribution from qq is larger than from gluon fusion. The total 
hadronic cross section in the MSSM is 22 times higher than in the SM (fa 175 pb). 

Out of all jets allowed by our cuts (11) a larger fraction of jets is radiated into the central 
part of the detector. For instance, the fraction of jets produced with |?73| < 2 compared to 
all jets allowed by the cuts is 93% in the MSSM versus 85% in the SM. Correspondingly, 
the pt spectrum is much softer than in the SM, yielding an enhanced rate for processes 
with jet transverse momenta below 850 GeV, e.g. by a factor of 10 for pt — 100 GeV, and 
rates similar to the SM above 850 GeV. 

4.3 Differential cross sections at the Tevatron 

The typical hadronic cross section for Higgs + jet in the SM expected at the Tevatron for the 
cuts Pt > 30 GeV and \rj3\ < 4.5 is around 0.1 pb for Higgs masses around 100 GeV, which 
is possibly not sufficient to be observable at the Tevatron. Therefore, for the Tevatron only 
the MSSM scenarios with low rrij^ and tan/3 not too small are of interest. Those scenarios 
exhibit a cross section enhanced by a factor of up to 30 compared to the SM [12]. This is 
due to the contribution of 6-quark initiated processes which become dominant because of 
the strongly enhanced Yukawa couphng of 6-quarks to the Higgs boson hP. 

Figure 9 shows results for the same low-m^ scenario as just described for the LHC in 
the previous paragraph. Very similar to the LHC case, we see a strongly enhanced total 
hadronic cross section with a softer pt spectrum and a larger fraction of jets radiated into 
the central part of the detector than in the SM. 
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5 Conclusions 



We have calculated pseudorapidity and transverse momentum distributions for the MSSM 
hP + high-pT jet production cross section at the LHC and the Tevatron. For scenarios with 
large m^, the loop-induced processes dominate the cross section, and superpartners can 
have a significant impact when they are not too heavy. For small m^, the Yukawa couplings 
of the 6-quarks are enhanced and hence the cross section is dominated by 6-quark induced 
tree-level parton reactions. The example investigated here, the m™^^(400) scenario, shows 
a strongly enhanced hadronic cross section compared to the SM, by a factor of more than 
20. Such a scenario predicts for both LHC and Tevatron a softer px spectrum, with a 
fraction of jets radiated into the central part of the detector larger than in the SM. 
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Figure 1: Typical quark and squark loop graphs for the process gg —>■ gh^ in leading order. 
Feynman graphs with opposite direction of charge flow are not depicted. 
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Figure 2: Typical Feynman graphs for the process ug uhP in leading order. For the 
scattering of the other quarks (c?, c, s) the graphs look similar. 
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Figure 3: Typical Feynman graphs for the process uu — gh^ in leading order. For the 
scattering of the other quarks {d, c, s) the Feynman graphs look similar. 
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Figure 4: Feynman graphs for the 6-quark processes in leading order. The graphs represent 
the amphtude for the process bg bh^ if the time axis points to the right and bb — > gh^ 
if it points down. 




Figure 5: LHC, no-mixing (700) scenario with itiao = 500 GeV, tan/3 = 35: hadronic cross 
section for Higgs + jet production differential in the jet's pseudorapidity 773 and transverse 
momentum pj- (lower left and right panel). Thick and thin lines indicate the MSSM and 
SM predictions respectively. Solid lines indicate the full result, dashed, dot-dashed and 
dotted lines the gg-, qg- and qq contribution respectively. In the upper panels the relative 
difference between the MSSM and SM result is displayed. 
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Figure 6: LHC, mh-max scenario with m^o = 400 GeV, tan/5 = 30 : differential hadronic 
cross sections for Higgs + jet production . See caption of Figure 5 for more details. 
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Figure 7: LHC, small alpha-eff scenario with m^o = 400 GeV, tan /3 = 30: differential 
hadronic cross sections for Higgs + jet production. See caption of Figure 5 for more 
details. 
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Figure 8: LHC, mh-max scenario with m^o = 110 GeV, tan/3 = 30: differential hadronic 
cross sections for Higgs + jet production. See caption of Figure 5 for more details. 
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Figure 9: Tevatron, mh-max scenario with m^o = 110 GeV, tan/5 = 30: differential 

hadronic cross sections for Higgs + jet production. See caption of Figure 5 for more 
details. 
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Figure 10: Relative and absolute difference between the MSSM and SM prediction for 
(Pa/drj^dpT at the LHC for the three benchmark scenarios as a function of pt and 773. The 
contour lines show the relative difference of the predictions in percent, while diamonds 
(♦), squares (■), triangles (a), circles (•), refer to an absolute difference in the range 0.1- 
0.5 fb/GeV, 0.5-1 flVGcV, 1-5 fb/GeV, 5-10 fb/GeV respectively. In the white area the 
difference is less than 0.1 fb/GeV. 
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